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Abstract There is compelling evidence for a strong electron-phonon interaction 
(EPI) in cuprate superconductors from the isotope effects on the su- 
percarrier mass, high resolution angle resolved photoemission spectro- 
scopies (ARPES), a number of optical and neutron-scattering measure- 
ments in accordance with our prediction of high-temperature super- 
conductivity in polaronic liquids. A number of observations point to 
the possibility that high- Tc cuprate superconductors may not be con- 
ventional Bardeen-Cooper-Schrieffer (BCS) superconductors, but rather 
derive from the Bose-Einstein condensation (BEC) of real-space pairs, 
which are mobile small bipolarons. Here I review the bipolaron theory 
of unconventional proximity effects, the symmetry and checkerboard 
modulations of the order parameter and quantum magneto-oscillations 
discovered recently in cuprates. 
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1. Polarons in high-temperature 
superconductors 

Many unconventional properties of cuprate superconductors may be 
attributed to the Bose-Einstein condensation (BEC) of real-space pairs, 
which are mobile small bipolarons [1, 2, 3, 4, 5, 6, 7, 8]. A possible funda- 
mental origin of such strong departure of the cuprates from conventional 
BCS behaviour is the unscreened Frohlich EPI providing the polaron 
level shift Ep of the order of 1 eV [4] , which is routinely neglected in the 
Hubbard U and t — J models. This huge interaction with c— axis polar- 
ized optical phonons is virtually unscreened at any doping of cuprates. 
Even acting alone it leads to large bipolarons in the continuum limit, if 
the ratio = e/eo of the high-frequency (electronic) and static dielectric 
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constants is small enough [3]. The large Frohlich bipolarons are further 
stabilised in going from 3D to 2D [9]. When acting together with the 
deformation potential and the Jahn- Teller EPIs the Frohlich EPI over- 
comes the inter-site Coulomb repulsion forming small bipolarons [5]. 
Hence, in order to build the adequate theory of high-temperature super- 
conductivity, the Coulomb repulsion and the unscreened long-range EPI 
together with the short-range one should be treated on an equal foot- 
ing with the short-range repulsive Hubbard U. When these interactions 
are strong compared with the kinetic energy of carriers, this "Coulomb- 
Frohlich" model predicts the ground state in the form of superlight small 
bipolarons [4, 10, 11]. 

Nowadays compelling evidence for a strong EPI has arrived from iso- 
tope effects [12, 13], more recent high resolution angle resolved pho- 
toemission spectroscopies (ARPES) [14], and a number of earlier optical 
[15, 16, 17, 18, 19], neutron-scattering [20] and recent inelastic scattering 
studies [21] of cuprates and related compounds. Whereas calculations 
based on the local spin-density approximation (LSDA) often predict neg- 
ligible EPI, the inclusion of Hubbard U in the LSDA + U calculations 
greatly enhances its strength [22]. 

A parameter-free estimate of the Fermi energy using the magnetic- 
field penetration depth found a very low value, ep < 100 nieV [23] 
clearly supporting the real-space (i.e individual) pairing in cuprate su- 
perconductors. There is strong experimental evidence for a gap in the 
normal-state electron density of states of cuprates [2] , which is known as 
the pseudogap. Experimentally measured pseudogaps of many cuprates 
are > 50 meV [24]. If following Ref. [25] one accepts that the pseudo- 
gap is about half of the pair binding energy. A, then the condition for 
real-space pairing, e^r < vrA, is well satisfied in most cuprates (typically 
the small bipolaron radius is rf, 0.2 — 0.4 nm). 

Also magnetotransport and thermal magnetotransport data strongly 
support preformed bosons in cuprates. In particular, many high-magnetic- 
field studies revealed a non-BCS upward curvature of the upper criti- 
cal field Hc2{T) (see [26] for a review of experimental data), in accor- 
dance with the theoretical prediction for the Bose-Einstein condensa- 
tion of charged bosons in the magnetic field [27]. The Lorenz number, 
L = e^K^/Ta differs significantly from the Sommerfeld value Lg = 7r^/3 
of the standard Fermi-liquid theory, if carriers are double-charged bosons 
[28]. Here Ke, and a are electron thermal and electrical conductivities, 
respectively. Ref. [28] predicted a rather low Lorenz number for bipo- 
larons, L « 0.15Le, due to the double elementary charge of bipolarons, 
and also due to their nearly classical distribution function above Tc- Di- 
rect measurements of the Lorenz number using the thermal Hall effect 
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[29] produced the value of L just above Tc about the same as predicted 
by the bipolaron model, and its strong temperature dependence. This 
breakdown of the Wiedemann-Franz law is apparently caused by ex- 
cited single polarons coexisting with bipolarons in the thermal equilib- 
rium [30, 31]. Also unusual normal state diamagnetism uncovered by 
torque magnetometery has been convincingly explained as the normal 
state (Landau) diamagnetism of charged bosons [32]. 

However, despite clear evidence for the existence of polarons in cuprates, 
no consensus currently exists concerning the microscopic mechanism of 
high-temperature superconductivity. While a number of early (1990s) 
and more recent studies prove that the Mott-Hubbard insulator pro- 
motes doping- induced polaron formation [33] , some other works suggest 
that EPI does not only not help, but hinder the pairing instability. The 
controversy should be resolved experimentally. Here I argue that the 
giant (GPE) and nil (NPE) proximity effects provide another piece of 
evidence for bipolaronic BEC in cuprates [34]. The same bipolaronic 
scenario also explains the symmetry and checkerboard modulations of 
the order parameter, and recently observed magnetooscillations in the 
vortex state of cuprates. 

2. Unconventional proximity effects 

Several groups [35] reported that in the Josephson cuprate SNS junc- 
tions supercurrent can run through normal A'^-barriers with the thickness 
2L greatly exceeding the coherence length, when the barrier is made from 
a slightly doped non-superconducting cuprate (the so-called N' barrier) , 
Fig.l. 

Using the advanced molecular beam epitaxy, Bozovic et al. [35] proved 
that GPE is intrinsic, rather than caused by any extrinsic inhomogene- 
ity of the barrier. Resonant scattering of soft-x-ray radiation did not 
find any signs of intrinsic inhomogeneity (such as charge stripes, charge- 
density waves, etc.) either [36]. Hence GPE defies the conventional ex- 
planation, which predicts that the critical current should exponentially 
decay with the characteristic length of about the coherence length, ^ ^ 1 
nm in cuprates. Annealing the junctions at low temperatures in vacuum 
rendered the barrier insulating. Remarkably when the SN'S junction 
was converted into a supcrconductor-insulator-superconductor (SIS) de- 
vice no supercurrent was observed, even in devices with the thinnest 
(one unit cell thick) barriers [37] (nil proximity effect, NPE). 

Both GPE and NPE can be broadly understood as the bipolaronic 
Bose-condensate tunnelling into a cuprate semiconductor [34]. 
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Figure 1. SNS cuprate nanostructure 



To illustrate the point one can apply the Gross-Pitaevskii (GP)-type 
equation for the superconducting order parameter V'(r), generalized by 
us [38] for a charged Bose liquid (CBL), since, as discussed above, many 
observations including a small coherence length point to a possibility 
that cuprate superconductors may not be conventional BCS supercon- 
ductors, but rather derive from BEG of real-space pairs, such as mobile 
small bipolarons [4], 



E{-ihV + 2eA) - n + j dr'y(r-r')|V'(r')P 



V'(r) = 0. (1) 



Here -E'(K) is the center-of-mass pair dispersion and the Peierls substi- 
tution, K —ihV + 2eA is applied with the vector potential A(r). 

The integro-differential equation (1) is quite different from the Ginzburg 
Landau [39] and Gross-Pitaevskii [40] equations, describing the order 
parameter in the BGS and neutral superfluids, respectively. Here ^ is 
the chemical potential and V{r) accounts for the long-range Goulomb 
and a short-range composed boson-boson repulsions [5] . While the elec- 
tric field potential can be found from the corresponding Poisson-like 
equation [38] , a solution of two coupled nonlinear differential equations 
for the order parameter il^{r) and for the potential V{r) in the nanos- 
tructure, Fig.l, is a nontrivial mathematical problem. For more trans- 
parency we restrict our analysis in this section by a short-range poten- 
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Figure 2. BEC order parameter at the SN boundary for jl = 1.0,0.1,0.01 and 
$J 0.001 (upper curve). The chemical potential is found above the boson band-edge 
due to the boson-boson repulsion in cuprate superconductors and below the edge in 
cuprate semiconductors with low doping. 

tial, V{r) = z;|^/;(r)p, where a constant v accounts for the short-range 
repulsion. Then in the absence of the magnetic field Eq.(l) is reduced 
to the familiar GP equation [40] using the continuum (effective mass) 
approximation, i^'(K) = K"^ /2mc. In the tunnelling geometry of SN' S 
junctions, Figs. 1,2, it takes the form, 



in the superconducting region, Z < 0, Fig. 2. Here rric is the boson 
mass in the direction of tunnelling along Z {h = c = ks = ^ this sec- 
tion). Deep inside the superconductor the order parameter is a constant, 
|^/;(Z)|2 = Us and ^ = vus , where the condensate density Ug is about 
Us ~ x/2, if the temperature is well below Tc of the superconducting 
electrode. Here the in-plane lattice constant a and the unit cell volume 
are taken as unity, and x is the doping level as in La2-a;Sra;Cu04. 

The normal barrier at Z > is an underdoped cuprate above its tran- 
sition temperature, < T where the chemical potential fi lies below the 
bosonic band by some energy e. Fig. 2, found from J dEN{E)[exp{{E -\- 
e)/r)-l]~i = x'/2. Here N{E) is the bipolaron density of states (DOS), 



1 d^^p{Z) 
2mc dZ"^ 



[HV(^)P-^]V(^) 



(2) 
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and x' < X is the doping level of the barrier. In-plane bipolarons are 
quasi-two dimensional repulsive bosons, propagating along the Cu02 
planes with the effective mass m several orders of magnitude smaller 
than their out-of-plane mass, nic 3> m, [5]. Using bipolaron band dis- 
persion, E(K) = K'^ /2m + 2tc[l — cos{K±d)], the density of states is 
found as N{E) = (m/27r2) arccos(l - E/2tc) for < E < At^ and 
N{E) = m/27r for 4tc < E. Here K and K± are the in-plane and out- 
of-plane center-of-mass momenta, respectively, tc = l/2mcd^, and d is 
the inter-plane distance. As a result one obtains 

e(r) ^ -rin(l-e-^"/^), (3) 

which is exponentially small at < T <C Tq turning into zero at T = T^, 
where « Tq/ ln(To/2tc), and Tq = vrx'/m 3> 3> tc- 

It is important to note that e{T) remains also small at T'^{2D) ^ T <C 
To in the purely two-dimensional repulsive Bose-gas [41]. While in two 
dimensions Bose condensation does not occur in either the ideal or the 
interacting system, there is a phase transition to a superfiuid state at 
T^{2D) = To/ln(l//o) ^ Tq, where /o <C 1 depends on the density 
of hard-core dilute bosons and their repulsion [42, 41]. The superfiuid 
transition takes place only if there is a residual repulsion between bosons, 
i.e. T^{2D) = for the ideal 2D Bose-gas. Actually T^{2D) gives a 
very good estimate for the exact Berezinski-Kosterlitz-Thouless (BKT) 
critical temperature in the dilute Bose gas, where the BKT contribution 
of vortices is important only very close to T^{2D) [42]. 

The GP equation in the barrier can be written as 

' ^'^(^)-[.|^(Z)p + e]^(Z). (4) 



2mr dZ^ 



Introducing the bulk coherence length, ^ = l/(2mcn^ju)^''^ and dimen- 
sionless f{z) = ip{Z)/nl^'^, p, = e/ugV, and z = Z/^ one obtains for a 
real f{z) 

if 2; < 0, and 

g = /' + A/, (6) 

if 2; > 0. These equations can be readily solved using first integrals of 
motion respecting the boundary conditions, /(— oo) = 1, and /(oo) = 0, 

£ = _(l/2 + /4/2_/2)i/2^ (7) 
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and 

= -{nf + f'm"\ (8) 



dz 

for z < and z > 0, respectively. The solution in the superconducting 
electrode is given by 



f{z) = tanh 



-2-/^.+o.5in^:;:i'+^!;;:+' 

2V2(l+/i)l/2 _1 



(9) 



It decays in the close vicinity of the barrier from 1 to /(O) = [2(l+/i)]^^/^ 
in the interval about the coherence length ^. On the other side of the 
boundary, z > 0, it is given by 

fu) = M!^^ (10) 

' sinh{z/iV2 + ln[2(/2(l + /i))V2 + (i + 4/i(i + /i))i/2]}- ^ ' 

Its profile is shown in Fig. 2. Remarkably, the order parameter pene- 
trates the normal layer up to the length Z* w which could be 
larger than ^ by many orders of magnitude, if ji is small. It is indeed 
the case, if the barrier layer is sufficiently doped. For example, taking 
x' = 0.1, c-axis rric = 2000me, in-plane m = lOrrie, a = 0.4 nm, and 
^ = 0.6 nm, yields Tq 140 K and > 50 at T = 25K. Hence 
the order parameter could penetrate the normal cuprate semiconduc- 
tor up to a hundred coherence lengths or even more as observed (GPE) 
[35]. If the thickness of the barrier L is small compared with Z*, and 
(/x)^/^ <C 1, the order parameter decays following the power law, rather 
than exponentially, 

m = (11) 

Hence, for L < Z*, the critical current should also decay following the 

power law rayher than exponentially. On the other hand, for the undoped 
barrier jl becomes larger than unity, jl cc \n.{mT /tix') oo for any finite 
temperature T when x' 0, and the current should exponentially decay 
with the characteristic length smaller that ^, which is experimentally 
observed as well (NPE) [37]. 

3. Quantum magneto-oscillations, d-wave 

symmetry and checkerboard modulations 

Until recently no convincing signatures of quantum magneto-oscillations 
have been found in the normal state of cuprate superconductors de- 
spite significant experimental efforts. There are no normal state oscil- 
lations even in high quality single crystals of overdoped cuprates like 
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Tl2Ba2Cu06, where conditions for de Haas-van Alphen (dHvA) and 
Shubnikov-de Haas (SdH) oscillations seem to be perfectly satisfied [43] 
and a large Fermi surface is identified in the angle-resolved photoe- 
mission spectra (ARPES) [44]. The recent observations of magneto- 
oscillations in kinetic [45, 46] and magnetic [47, 48] response functions 
of underdoped YBa2Cu306.5 and YBa2Cu408 are perhaps even more 
striking since many probes of underdoped cuprates including ARPES 
[49] clearly point to a non Fermi-liquid normal state. Their descrip- 
tion in the framework of the standard theory for a metal [50] has led 
to a very small Fermi-surface area of a few percent of the first Brillouin 
zone [45, 46, 47, 48], and to a low Fermi energy of only about the room 
temperature [47]. Clearly such oscillations are incompatible with the 
first-principle (LDA) band structures of cuprates, but might be compat- 
ible with a a low Fermi energy and non-adiabatic polaronic normal state 
of charge-transfer Mott insulators as discussed in section 1. 

Nevertheless one can raise a doubt concerning their normal state ori- 
gin. The magnetic length, A = {Trh/eB)^^"^ > 5 nm, remains larger 
than the zero-temperature in-plane coherence length, ^ < 2 nm, mea- 
sured independently, in any field reached in Ref. [45, 46, 47, 48]. Hence 
the magneto-oscillations are observed in the vortex (mixed) state well 
below the upper critical field, rather than in the normal state, as also 
confirmed by the negative sign of the Hall resistance [45]. It is well 
known, that in "YBCO" the Hall conductivities of vortexes and quasi- 
particles have opposite sign causing the sign change in the Hall effect 
in the mixed state [51]. Also there is a substantial magnetoresistance 
[46] , which is a signature of the flux flow regime rather than of the nor- 
mal state. Hence it would be rather implausible if such oscillations have 
a normal-state origin due to small electron Fermi surface pockets [48] 
with the characteristic wave-length of electrons larger than the widely 
accepted coherence length. 

Here I propose an alternative explanation of the magneto-oscillations 
[45, 46, 47, 48] as emerging from the quantum interference of the vortex 
lattice and the checkerboard or lattice modulations of the order param- 
eter observed by STM with atomic resolution [52]. The checkerboard 
effectively pins the vortex lattice, when the period of the latter, A is 
commensurate with the period of the checkerboard lattice, a. The con- 
dition A = Na , where N is a large integer, yields l/B^^"^ periodicity 
of the response functions, rather than 1/B periodicity of conventional 
normal state magneto-oscillations. 

The integro-differential equation (1) in the continuum approximation, 
£'(K) = h'^K^/2m, with the long-range Coulomb repulsion between dou- 
ble charged bosons, V{r) = K(r) = Ae'^/eor, describes a single vortex 
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Figure 3. The order parameter profile f{p) = %p{r) /ri's^'^ of a single vortex in CBL 
[38] (symbols). Inset: CBL vortex (a) [38, 54] compared with the Abrikosov vortex 
(b) [53] {herep=[x^ +y^f'^). 



with a charged core, Fig. 3, when the magnetic field, B is apphed. The 
coherence length in this case is roughly the same as the screening radius, 
^ = {h/2^/'^mujpY/'^ . Here ujp = {IQ-KUge^ / CQm)^/'^ is the CBL plasma 
frequency, eg the static dielectric constant of the host lattice, m is the 
(in-plane) boson mass, and is the average condensate density. The 
chemical potential is zero, = 0, if one takes into account the Coulomb 
interaction alone due to a neutralizing homogeneous charge background. 
Each vortex carries one flux quantum, (j)Q = nh/e , but it has an unusual 
core, Fig.Sa, Ref. [38], due to a local charge redistribution caused by 
the magnetic field, different from the conventional vortex [53], Fig. 3b. 
Remarkably, the coherence length turns out very small, ^ ~ 0.5nm with 
the material parameters typical for underdoped cuprates, m = lOmg, 
Hs = lO^^cm"^ and eo = 100. 

The coherence length ^ is so small at low temperatures, that the 
distance between two vortices remains large compared with the vortex 
size, A » ^, in any laboratory field reached so far [45, 46, 47, 48]. 
It allows us to write down the vortex-lattice order parameter, 'il>{r) = 
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ipvliT^), as 



n 



1/2 



(12) 



where ^(r) = 1 — f{p), and = X{nx,ny} with nx,y = 0, ±1,±2, ... (if, 
for simphcity, we take the square vortex lattice [55]). The function (f){p) 
is Hnear well inside the core, (p^p) « 1 — 1.52p/^ (p <C ^), and it has a 
small negative tail, (/>(p) « —4^^/ outside the core when p ^ ^, Fig.3 
[38]. 

In the continuum approximation with the Coulomb interaction alone 

the magnetization of CBL follows the standard logarithmic law, M{B) cx 
In 1/5 without any oscillations since the magnetic field profile is the same 
as in the conventional vortex lattice [54]. However, more often than not 
the center-of-mass Bloch band of preformed pairs, £^(K), has its minima 
at some finite wave vectors K = G of their center-of-mass Brillouin zone 
[4, 5]. Near the minima the GP equation (1) is written as 



(-iW -hG + 2eAf 

2m** 



V'(r)+ 



J dr'V{r - r')|V'(r')PV'(r) = 0, 



(13) 



with the solution ■0(r) = 'i/'G(r) = e^^'^ijjyi{r), if the interaction is the 
long-range Coulomb one, V{r) = lc(r). 

In particular, a nearest-neighbor (nn) approximation for the hopping 
of intersite bipolarons between oxygen p-orbitals on the Cu02 2D lattice 
yields four generate states ipa with Gj = {±27r/ao, ±27r/ao} , where Oq 
is the lattice period [10]. Their positions in the Brillouin zone move 
towards F point beyond the nn approximation. The true ground state is 
a superposition of four degenerate states, respecting time-reversal and 
parity symmetries [56], 



i/j{r) = An]^'^ [cos(7ra;/a) ± cos(7ry/a)] i/^vli'"^)- 



(14) 



Here we use the reference frame with x and y axes along the nodal di- 
rections and a = 2^^/^ao. Two "plus/minus" coherent states, Eq.(14), 
are physically identical since they are related via a translation trans- 
formation, y ^ y + a. Normalizing the order parameter by its average 
value (V'(r)^) = Ug and using (^/A)^ <C 1 as a small parameter yield the 
following "minus" state amplitude, A a 1 - iVE^o 2['^i(2^/^vr/a) + 

(f>2{2^^'^ir/a)]5n^R/2 + [(t>ii2TT / a) + (f>2{2iT / a)]6n,R for the square vortex lat- 
tice [55] with the reciprocal vectors g = {27r /X){nx,ny}. Here 6n,R 
is the Kroneker symbol, R = A/a is the ratio of the vortex lattice 
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Figure 4- The checkerboard d-wave order parameter of CBL [56] on the square 
lattice in zero magnetic field (coordinates x,y are measured in units of a). 



period to the checkerboard period (n = 0,1,2,...), = BS/cpo is 
the number of flux quanta in the area S of the sample, and (pkiQ) = 
(27r/5) dppJo{pq)(p^ (p) is the Fourier transform of A:'s power of (t){p), 
where Jq{x) is the zero-order Bessel function. 

The order parameter ip{r), Eq.(14) has the d-wave symmetry changing 
sign in real space, when the lattice is rotated by 7r/2. This symmetry is 
due to the pair center-of-mass energy dispersion with the four minima at 
K 7^ 0, rather than due to a specific symmetry of the pairing potential. 
It also reveals itself as a checkerboard modulation of the carrier density 
with two-dimensional patterns in zero magnetic field. Fig. 4, as predicted 
by us [56] prior to their observations [52]. Solving the Bogoliubov-de 
Gennes equations with the order parameter, Eq.(14), yields the real- 
space checkerboard modulations of the single-particle density of states 
[56], similar to those observed by STM in cuprate superconductors. 

Now we take into account that the interaction between composed 
pairs includes a short-range repulsion along with the long-range Coulomb 
one, V{r) = V^(r) + v5{r) [5]. At sufficiently low carrier density the 
short-range repulsion can be treated as a perturbation to the ground 
state, Eq.(14). Importantly the short-range repulsion energy of CBL, 
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U = {v/2){ip{r)^), has a part, AU, oscillating with the magnetic field as 
AU °° 

iV^ [Ai5^^R/2+A26n,R + A 



^0 



where ?7o = vn1/2 is the hard-core energy of a homogeneous CBL, and 
the amplitudes are proportional to the Fourier transforms of 0(p) as 

Ai = 15^1 (^2V27r/a) - 45^2 {2}''^TT/a^ + 24^3 (^2VV/a) - 6^4 (2^^^ 
+ 8(^1 (^lO^V/o) - 1202 (10^/2^/a) + 8<^3 (lO^/^^/a) - 2(^4 (lO^/V/a) 



= -(23/2)^1 (27r/a) + (57/2)(^2 (27r/a) - 16^3 (27r/a) + 4^4 (27r/a) 
- 12^1 (|2^/V/a) + 9^2 (2=^/V/a) - 6^3 (2=^/V/a) + 8^4 (2=^/V/a) 

(17) 

A3 = -01 (47r/a) + (3/2)02 (47r/a) - 03 (4vr/a) + (1/4)04 (4^/a) . (18) 

Fluctuations of the pulsed magnetic field and unavoidable disorder in 
cuprates induce some random distribution of the vortex-lattice period, 
A. Hence one has to average AU over R with the Gaussian distribution, 
G{R) = exp[-(i? - i?)V72]/77ri/2 around an average R with the width 
^ <^ R. Then using the Poisson summation formula yields 

^ = iV y Aie-'^'^V/i6 ,os(^fci?) 

+^26-"'^='^'/^ cos(27r A-i?) + ^36-'^'*^'^' cos(47rA:i?) . (19) 

The oscillating correction to the magnetic susceptibility, Ax{B) = 
—d^Cl/dB^, is strongly enhanced due to high oscillating frequencies in 
Eq.(19). Since the superfluid has no entropy we can use AU as the 
quantum correction to the thermodynamic potential even at finite 
temperatures below Tc{B). Differentiating twice the first harmonic {k = 
1) of the first lesser damped term in Eq.(19) we obtain 

Ax{B) ^ Xoe-''^°/''^ (§) ' cos(Bo/5)V2, (20) 

where xo = UoSAie^a^/Air^h'^ is a temperature-dependent amplitude, 
proportional to the condensate density squared, Bq = n^h/ea^ = Sir^h/eaQ 
is a characteristic magnetic field, which is approximately 1.1 • 10^ Tesla 
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Figure 5. Quantum corrections to the vortex- lattice susceptibility versus 1/-B, 
Eq.(20) (solid line, Bo = 1.000- lO'^ Tesla, 5 = 0.06) compared with oscillating suscep- 
tibility of YBa2Cu306.5 (symbols) and with the conventional normal state oscillations 
(dashed line) [48] at T = 0.4 K. 



for ao ~ 0.38 nm, and 7 is replaced by 7 = 6R with the relative distribu- 
tion width 5. Assuming that C ^ a, so that the amplitude Ai is roughly 
a'^/S, the quantum correction A^, Eq.(20), is of the order of wx'^/B'^, 
where x is the density of holes per unit cell. It is smaller than the conven- 
tional normal state (de Haas-van Alphen) correction, AxdHvA ~ t^/B"^ 
[50], for a comparable Fermi-energy scale /i = wx, since x <C 1 in the 
underdoped cuprates. 

Different from normal state dHvA oscillations, which are periodic 
versus 1/-B, the vortex-lattice oscillations, Eq.(20) are periodic versus 
l/B^/"^. They are quasi-periodic versus \/B with a field-dependent fre- 
quency F = Bq{B / BqY^'^ /2-K , which is strongly reduced relative to the 
conventional-metal frequency (~ Bq/I-k) since B <C Bq, as observed in 
the experiments [45, 46, 47, 48]. The quantum correction to the sus- 
ceptibility, Eq.(20) fits well the oscillations in YBa2Cu306.5 [48], Fig.5. 
The oscillations amplitudes, proportional to exp(— 5^i?o/16-B) decay 
with increasing temperature since the randomness of the vortex lattice, 
5, increases, and the Bose-condensate evaporates. 
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4. Summary 

A possibility of real-space pairing, proposed originally by Ogg [57] 
and later on by Schafroth and Blatt and Butler [58] has been the sub- 
ject of many discussions as opposed to the Cooper pairing, particularly 
heated over the last 20 years after the discovery of high temperature su- 
perconductivity in cuprates. Our extension of the BCS theory towards 
the strong interaction between electrons and ion vibrations [5] proved 
that BCS and Ogg-Schafroth pictures arc two extreme limits of the same 
problem. Cuprates are characterized by poor screening of high-frequency 
optical phonons, which allows the long-range Prohlich EPI bound holes 
in siipcrlight small bipolarons [4], which are several orders of magnitude 
lighter than small Holstein (bi)polarons. 

The bipolaron theory accounts for GPE and NPE in slightly doped 
semiconducting and undoped insulating cuprates, respectively. It pre- 
dicts the occurrence of a new length scale, H/ y^2mce{T), which turns 
out much larger than the zero-temperature coherence length in a wide 
temperature range above the transition point of the normal barrier, if 
bosons arc almost two-dimensional (2D). The physical reason, why the 
quasi-2D bosons display a large normal-state coherence length, whereas 
3D Bose-systems (or any-D Fermi-systems) at the same values of pa- 
rameters do not, originates in the large DOS near the band edge of two- 
dimensional bosons compared with 3D DOS. Since DOS is large, the 
chemical potential is pinned near the edge with the magnitude, e(T), 
which is exponentially small 

Here I also propose that the magneto-oscillations in underdoped cuprate 
superconductors result from the quantum interference of the vortex lat- 
tice and the lattice modulations of the d-wave order parameter, Fig.4, 
which play the role of a periodic pinning grid. Our expression (20) de- 
scribes the oscillations as well as the standard Lifshitz-Koscvich formula 
of dHvA and SdH effects [45, 46, 47, 48]. The difference of these two 
dependencies could be resolved in ultrahigh magnetic fields as shown 
in Fig. 5. While our theory utilizes GP-type equation for composed 
charged bosons [38], the quantum interference of vortex and crystal lat- 
tice modulations of the order parameter is quite universal extending well 
beyond Eq.(l) independent of a particular pairing mechanism. It can 
also take place in the standard BCS superconductivity at -B < Hc2, 
but hardly be observed because of much lower value of Hc2 in con- 
ventional superconductors resulting in a very small damping factor, 
oc exp{-5'^Bo/16B) <^ 1. 

I appreciate valuable discussions with A. F. Andreev, I. Bozovic, 
A. Bussmann-Holder, J. T. Devreese, L. P. Gor'kov, V. V. Kabanov, 
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and support of this work by EPSRC (UK) (grant Nos. EP/D035589, 
EP/C518365). 
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